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ABSTRACT 

The space velocity of a stellar black hole encodes the history of its formation 
and evolution. Here we measure the 3-dimensional motion of the microquasar 
GRS 1915+105, using a decade of astrometry with the NRAO Very Long Baseline 
Array, together with the published radial velocity. The velocity in the Galactic 
Plane deviates from circular rotation by 53-80±8 km s _1 , where the range covers 
any specific distance from 6-12 kpc. Perpendicular to the plane, the velocity is 
only 10±4 km s _1 . The peculiar velocity is minimized at a distance 9-10 kpc, 
and is then nearly in the radial direction towards the Galactic Center. We discuss 



1 On leave from Service d'Astrophysique, CEA-Saclay, 91191 Gif-sur-Yvette, France. 
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mechanisms for the origin of the peculiar velocity, and conclude that it is most 
likely a consequence of Galactic velocity diffusion on this old binary, rather than 
the result of a supernova kick during the formation of the 14 M black hole. 
Finally, a brief comparison is made with 4 other BH binaries whose kinematics 
are well determined. 

Subject headings: astrometry — stars: kinematics — stars: individual (GRS 1915+105) 
- techniques: high angular resolution — techniques: interferometric 



Introduction 



GRS 1915+105 is a famous microquasar with superluminal ejecta ( Mirabel fc Rodriguez! 
1994), and spectac ular X-ray variability from an unstable accretion disk (jGreiner et al.lll996l : 
Belloni et al.ll2000i ). For nearly a decade, its natu re was highly obscur ed at a distance D of 
~10 kpc. Finally, IR spectroscopy with the VLT (jGreiner et al.ll200ll ) was used to measure 
the heliocentric radial velocity of —3+10 km s _1 for the center of mass of the binary, a 
mass of 14+4 M for the compact o bject, with a donor m ass of ~1.2±0.2 M in an orbit 
of 33.5 days period. VLBA imaging (IDhawan et al.ll2000al ). in addition to the ejecta during 
outbursts, revealed an AU-sized radio jet closely coupled to the accretion disk activity. Under 
the conditions imaged here, the so-called low-hard state, the flux of the AU-size jet follows 



the va riation in X-rays with a time-lag of a few minutes, (IMirabel et al.lll998l ; IDhawan et al. 
2000al ). The jet is hence located within 1 AU (0.1 mas) of the black hole. 



2. Observations and Results 



2.1. Observations 



Our observations spanned ~10 years, sampled at random times dictated by the flare 
activity of the source. Usually the goal was to image the AU-size jet and larger scale 
superluminal ejecta in coordination with X-ray instruments. The moving ejecta can confuse 
the position determination. Indeed there were some occasions just after a major flare when 
the central nucleus was very weak or absent, and the moving ejecta dominated the flux 
density of the source. We have hence selected for this astrometric study only those images 
which were dominated by the central source of flux density 10 to 50 mJy, co rresponding to 



(2000a 


); 


Fuchs et al. 


(2003 


); 


Ribo et al. 


(2004 
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Standard VLB A phase- referencing was used, (IBeasley h Conwaylll995l ) with a total cy- 
cle time of 3 minutes, of which 1 minute was on the primary calibrator J1924+1540, at a 
separation of 5.3° from the target. This technique provides the position of the black hole 
relative to the reference source, whose position has been pre-determined in the ICRF frame. 
Most of the observations were at a frequency of 8.4 GHz (3.6cm wavelength). During some 
observations, depending on the radio spectra, (if available from monitoring with other in- 
struments) we observed at 3.6 and 2 cm. These were alternated every 20min, at 256Mbit s _1 
recording. Standard post-processing corrections were applied for ionospheric refraction us- 
ing a GPS-based global model. A-priori amplitude calibration was done from measured 
system temperatures and antenna gains, followed by self-calibration on the reference source, 
application of the refined calibration to all sources including the target, and imaging. 

We fitted a single elliptical Gaussian to each image to extract the position. The formal 
errors of the fit are far exceeded by tropospheric variations, resulting in position errors 
~1 mas rms at each epoch. The earliest observations have larger error due to the poor 
choice of reference sources, J1922+1530 (4.8 away, but quite resolved) and J1925+2106, 
10.5° away from the target. We estimated the errors, including systematic errors peculiar to 
the weather on that day, from the scatter in image position over intervals of ~20 min. In 
later observations, a second calibrator, J1922+0840 (2.9 ° away), was observed as a quality 
check and processed in identical fashion to the target. Positions at the two wavelengths on 
a given observation date agreed, and were averaged together for that date. On one occasion 
observations were possible on successive days, positions from which were averaged together. 



2.2. Results 

The measurements Q are summarized in Table-[H and Figure-[I1 showing the coordinates 
of the black hole versus time. We fit the following position and proper motion at MJD 51545 
(01 January 2000): 

RA (J2000) 19 h 15 m ll s .54902 ±0.7 mas Dec(J2000) +10°56'44" .7478 ±0.9 mas 
/Xra Cos(Dec) = -2.86 ±0.07 mas yr" 1 // Dec = -6.20 ±0.09 mas yr" 1 

1 Strictly, the measured VLBA positions are Geocentric. In the table, plot & fit, we report measured 
quantities, rather than correcting them for parallax with an uncertain distance. In the fit for proper motion, 
the error is reduced by the random observation date, and the long timespan. The parallax is 0.1 mas at 
_D=10 kpc. The changes in proper motion if we correct for parallax are ~0.01 mas yr _1 , far below our 
measurement errors. Thus Geocentric positions can safely be taken as Heliocentric. 
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The corresponding Heliocentric motions in Galactic coordinates are: 



-6.82 ±0.10 mas yr 



-i 



-0.35 ±0.07 mas yr" 1 



Using our measured proper motion and the known V ra di a i of — 3±1 km s" 1 , we transform 
these measured Heliocentric quantities with the standard method of I Johnson fc Soderblom 
(119871 ) updated for J2000, and the Solar peculiar motion mea sured by Hipparcos, [U , V , 
V = 10.0, 5.25, 7.17 km s" 1 ] from iDehnen fc Binnevl (119981 ). We use the convention U ± 
toward I = 0°; V + toward I = 90°; W + toward the North Galactic Pole. 

We compare these measured velocities with those expected from circular rotation at 
220 km s~ l at the locat i on of the black hole. For the distance, 11 to 12 kpc have been 



proposed, (iFenderl 120031 ; iDhawan et a 



suggested by others, e.g. Kaiser et al 



2000 



a), but since considerably smaller values are 



(120041 ). we let the distance be a free parameter in 



Figure [2] and Table [2j We then evaluate the velocity discrepancy in local and Galactocentric 
frames, and note that: 



The measured W component is small, so the black hole velocity lies primarily in the 
Galactic Plane. 

At any distance there is a significant component of measured velocity directed inward 
to the Galactic Center, whereas none is expected for purely circular motion. 

The extremum (negative) proper motion from Galactic rotation along this line of sight 
occurs at D— 9 kpc, and is —5.85 mas yr"" 1 , while —6.82 mas yr"" 1 is measured. 

Figure [3] shows that the total peculiar velocity of the black hole (in either LSR or GC 
frame) is minimized at D— 9-10 kpc, and is then 53 km s -1 in the radial direction 
towards the Galactic Center. 

The measured heliocentric vel ocity matches t hat e xpected from circular rotation at 
~12 kpc, as already noted by iGreiner et al.l (120011 ). However, at this distance, the 
total peculiar velocity is 83 km s _1 . 

At D= 6.5 kpc, both the Galactocentric components mismatch circular rotation by 
over 45 km s _1 each. 



Hence, the distance favoured by velocity matching is 9-10 kpc, but this is not a very 
strong argument. A distance determined by trigonometric parallax, while possible with the 
VLBA, remains very challenging for this object. 



- 5 - 



2.3. The Galactocentric Orbit 



We next compute the Galactocentric orbit of the blac k hole assuming a smoo th Galactic 
potential composed of disk, spheroid and halo, following IJohnston et al.l (119951 ). Figure 0] 
shows the orbit assuming a distance of 12 kpc. The path for the past 1 Gyr is the blue 
(dashed) line, with the last 250 Myr in red (solid line). The orbit deviates from circularity 
with an eccentricity of 0.285, and reaches a maximum height above the Galactic Plane of 
only 230 pc, typical of the thin disk population. However, GRS 1915+105 does not match 
the velocity of thin disk objects at any distance. Detailed computation of the orbit including 
the effects of spiral arms is beyond the scope of this paper. 



3. Discussion 

We discuss mechanisms by which the black hole could have acquired its anomalous 
velocity. To summarize, though the velocity is significant, its component out of the Galactic 
plane is low, so an impulsive origin cannot be proved. The age and uncertain distance of 
the system disallow a strong conclusion about the origin of the velocity. However, our data 
helps to rank the probability of three mechanisms. 



3.1. Stellar Diffusion 



An individual star in the disk can acquire significant changes in angular momentum 
from non-axisym metric forces such as spiral arms, leadin g to an increase of its non-circular 
motion with age (IWeilenl 119771 ; ISellwood fc Binneyl 120021 ) . If we assume the binary is as old 
as the donor, a red giant of 1.2 M Q aged ~4 Gyr, then it has endured about 16 Galactic 
orbits. Velocity imparted by spiral arm passages conforms naturally to the measured low 
velocity perpendicular to the Galactic plane, and has had sufficient time to operate, given 
the age of the system. In comparison, the non-circular velocity of progres sively old e r disk 



stars in the so lar neighborhood i s 30-50 km s _1 , as discussed for example in iMignardl (120001 . 
Table 6), and lAlcobe fc Cubarsil (120051 . population components S3 and A 3 ). At a distance 
D ~9-10 kpc, the above mechanism suffices to account for the peculiar velocity, and must 
have operated to some extent. It may have difficulty in producing a velocity perturbation of 
70-80 km s _1 , as would be the case if the source were at 6 or at 12 kpc. 
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3.2. Natal Kicks 

We note that even at 9 kpc, the black hole possesses the momentum of a pulsar of 
1.4 M Q with velocity of 530 km s _1 . Hence in the following we discuss the possible origin of 
the peculiar velocity in a natal supernova event, assuming one occurred. We examine this 
assumption in a later section. 



3.2.1. Blaauw Kick 



Mass unbound suddenly and symmetrically by the supernova explosion forces the binary 
to recoil with a momentum opposite to that of the unbound mass at the time of ejection. 
This spherical mass loss resul ts in a kick that lies i n the plane of the binary orbit. We 
estimate, using equation (1) of iNelemans et al.l (119991 ). that mass unbound from the binary 
in a SN explosio n impa rted at most 23-44 km s _1 , for allowable mass ratios as measured 
by iGreiner et al.l (120011 ). Assuming the binary plane was randomly oriented, the kick has 
~25% probability of lying in the galactic plane within ±15°, which is the ratio of measured 
peculiar velocities perpendicular and parallel to the plane. 



3.2.2. Asymmetric Kick: 



If the binary lies at D— 12 kpc or 6 kpc, an additional kick of ~30 km s _1 is re- 
quired. It is generally accepted that neutron stars form in asymmetric SN, and receive 
substantial birth ki cks, which are seen in the velocities of young single pulsars, ~400km s~ ] 



( Hobbs et al. 2005h . Neutron stars in long period binaries may receive smaller kicks of 



~50 km s _1 (IPfahl et al.ll2002l ). which is a case more closely related to the present binary 
with orbital period 33 days. This kick due to any asymmetries in the explosion has no pre- 
ferred direction relative to the binary orbit plane, unless it is assumed that the progenitor's 
spin is parallel to the orbital axis. Assuming random orientation, an asymmetric supernova 
has roughly the same 25% chance (as above), of producing a kick that lies in the Galactic 
plane. 



3.3. Comparison with the kinematics of other black hole binaries 



The question of whether or not black holes receive a natal kick similar to neutron stars 
is an open one. On the theoretical side, the formation of black holes of over 10 M Q challenges 
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stellar evolution models. iFryer fc Kalogeral (120011 ) proposed that black holes of that mass 
should form by direct collaps e without SN explos i on, or after a SN with energy insufficient to 
unbind the stellar envelope. iBelczynski fc Bulikl (120021 ) also conclude that standard models 
would work only if wind mass loss for very massive stars was half as much as commonly 
assumed, and black holes of >10 M would form by direct collapse. Hence no kick velocity 
is expected in the formation of the most massive BH. There is some observational evidence 
that energetic SN may not occur in the formation of stellar black holes of relatively high 
mass. Cyg X-l is a 10 M Q B H with a peculiar velocity o f <10 km s _1 , and may be an 
example of such a 'dark birth' (IMirabel fc Rodriguesll2003al ). 



Less massive black holes can form through partial fallback onto the collapsing object, 
after a SN explosion. We briefly evaluate the evidence for runaway black holes which have 
well measured, 3-dimensional velocities. The most likely mechanism for their high velocity 
is an asymmetric kick from a SN. Independent confirmation comes from rather difficult mea- 
surement of the supernova's nucleos ynthesis products in the companio n star. Such evidence 
exists for neutron stars as well, e.g., IGonzalez Hernandez et al.l (120051 ). 



The runaway b inary GRO 1655—40 contains a dynamically con firmed BH of 6.3 

and has an anomalous velocity v= 120 km s 1 ( IMirabel et al.ll2002l ). 



(jGreene et al 



Israelian et al.l (119991 ) found super-solar abundances of several a-elements in the atmo- 
sphere of the compa nion star, and concluded that SN ejecta had enriched the companion. 
Foellmi et al.l (120071 ), using different techniques, found little evidence for enrichment except 



in O. More recently, IGonzalez Hernandez et al.l (120071 ) . using VLT/UVES spectra, support 



the abundances derived by 



Israelian et al 



(J1999|). 



XTE 1118+48 is another confirmed 7-8 M^ BH binary (Gelino et al.ll2006l ). orbiting 



in the Galactic halo with v= 145 km s 1 . IMirabel et al.l (120011) proposed that it was born 



in a globular cluster ~7 Gyr ago. However, more recently IGonzalez Hernandez et al.l (120061 ) 
claim to have found fossil evidence of a supernova in the metal enrichment of the small mass 
donor, and propose instead that the progenitor was born in the Galactic disk and the binary 
shot out to the halo by the SN explosion. 



LS 5039 has a runaway velocity of 126 km s" 1 , (IRibo et al.l 120021 ). and has clearly 
been kicked, but its mass estimate is closer to the ne utron star range th an the objects 
discussed above. It may contain a black hole of 3.7M (jCasares et al.ll2005l ). if the binary 
is synchronized at periastron. Without the assumption of pseudo-synchroniz ation, the mass 



might be as low as 1.4M . Regarding the direct detection of SN byproducts, iMcSwain et al. 



( 120041 ) find N enrichment, but they cautio n that it could result from internal mixing due to 
rapid rotation, a scenario also favoured by lCasares et al.l (120051 ). 
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Finally, GRS 1915+105 does not fit easily in either the high or low speed category. 
The velocity is significantly higher than Cyg X-l. However, we cannot conclude that 
GRS 1915+105 is a runaway, because its low velocity out of the plane, and the operation 
of velocity diffusion in the plane, imply that any natal kick (symmetric and/or asymmetric) 
was pro bably smaller t h an ~ 30 km -1 . We hope our data will permit an analysis such as 
done by IWillems et al.l (120051 ) for GRO J1655— 40, who used all available information to 
reconstruct the evolution back to the time of core collapse and compact object formation. 
The observations, even for the small sample so far, imply a variety of formation scenarios 
for black holes. 



4. Conclusions 

1. The most massive known stellar black hole moves at (53-80) ±8 km s _1 with respect 
to its regional standard of rest, for any specific distance in the range 6-12 kpc. 

2. At the distance of 9-10 kpc, the peculiar velocity is a minimum. If GRS 1915 is 3- 
5 Gyr old, it could have acquired this velocity by dynamic diffusion, i.e., non-circular 
motion imparted by perturbations such as spiral arms over time. In this case the black 
hole in GRS 1915+105 would have formed without a natal kick, like the 10 M black 
hole in Cyg X-l. 

3. If the distance is smaller (6 kpc) or larger (12 kpc), then some additional ~30 km s _1 
is required from another mechanism such as a natal supernova explosion, and the kick 
is constrained to lie within ±15° of the Galactic plane. 
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Table I. ICRF Position of Nucleus of GRS 1915+105. 



MJD 




RA 


CRA 


Dec 


CDec 




Comment 


50227 


0. 


.54975 


1, 


.8 





.7685 


1 


.8 


1996 


May 24 / 8.4 GHz 


50229 


0. 


.54970 


1 


.8 





.7706 


1 


.8 


1996 


May 26 / 8.4 GHz 


50299 


o. 


.54965 


1, 


.5 





.7670 


1 


.5 


1996 


Aug 03 / 8.4 GHz 


50583 


0. 


.54955 


1, 


,5 





.7634 


1 


.5 


1997 


May 15 / 8.4, 15 GHz 


50744 


0. 


.54943 


1, 


.5 





.7621 


1 


.5 


1997 


Oct 23 / 15 GHz 


50752 





.54945 


1, 


.5 





.7629 


1 


.5 


1997 


Oct 31 / 8.4 GHz 


50913 





.54939 


1, 


.0 





.7589 


1 


.0 


1998 


Apr 11 / 15 GHz 


50935 





.54937 


1. 


.0 





.7582 


1 


.0 


1998 


May 02/15 GHz 


51331 





.54906 


1, 


.0 





.7522 


1 


.0 


1999 


Jun 02 / 8.4, 15 GHz 


51338 





.54905 


1. 


.0 





.7523 


1 


.0 


1999 


Jun 09 / 8.4 GHz 


52107 


0. 


.54870 


0. 


.8 





.7384 





.8 


2001 


Jul 17,18 / 8.4,15 GHz 


52722 





.54846 


0. 


.5 





.7269 


0. 


.7 


2003 


Mar 24 / 8.4, 15 GHz 


52731 


0. 


.54843 


0. 


.7 





.7265 





.9 


2003 


Apr 02 / 8.4, 15 GHz 


52748 





.54840 


0. 


.5 





.7264 





.7 


2003 


Apr 19 / 8.4, 15 GHz 


53794 


0. 


.54778 


0. 


.7 





.7093 





.8 


2006 


Feb 28 / 8.4 GHz 


53800 


0. 


.54777 


0. 


.7 





.7095 





.8 


2006 


Mar 06 / 8.4 GHz 



Note. — - 

RA(J2000) offset in seconds from 19 h 15 m ll s . 
ora in milliarcseconds, i.e. multiplied by Cos(Dec) 
Dec(J2000) offset from 10°56'44" 
a"Dec in milliarcseconds. 
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Table 2. Velocity Components of GRS 1915+105 vs. Distance. 



D 


U 






V 




W 


aw 


UC 


VC 








VR G 


VC G 


Vp 


5 


122. 


.9 


7.5 


-110.5 


7.6 


-1.1 


1.8 


127.8 


-40. 


.9 


49.9 


-116.2 


-36.4 


160.5 


69.8 


6 


145 


.9 


7.6 


-133.2 


7.7 


-2.8 


2.1 


155.3 


-64. 


.2 


52.7 


-152.2 


-42.1 


164.5 


69.7 


7 


168 


.9 


7.7 


-155.9 


7.7 


-4.4 


2.5 


178.6 


-91. 


.6 


49.6 


-188.0 


-46.6 


174.5 


65.2 


8 


191 


.9 


7.7 


-178.6 


7.8 


-6.1 


2.8 


196.3 


-120. 


.7 


41.3 


-221.1 


-49.7 


189.9 


58.4 


9 


214 


.9 


7.8 


-201.4 


7.9 


-7.7 


3.2 


208.3 


-149 


.2 


29.4 


-249.6 


-51.5 


209.5 


53.2 


10 


237 


.9 


7.9 


-224.1 


8.0 


-9.4 


3.5 


215.4 


-175 


.4 


15.8 


-273.1 


-52.3 


232.1 


54.5 


11 


260 


.9 


8.1 


-246.8 


8.1 


-11.1 


3.9 


218.9 


-198 


.4 


1.9 


-291.8 


-52.3 


257.0 


65.0 


12 


283 


.9 


8.2 


-269.5 


8.2 


-12.7 


4.3 


220.0 


-218 


.2 


-11.5 


-306.4 


-51.8 


283.5 


82.9 


13 


306 


.9 


8.3 


-292.2 


8.3 


-14.4 


4.6 


219.5 


-235 


.0 


-23.8 


-317.9 


-51.1 


311.1 


105.4 
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Note. — Column headings are: 
1: D = Distance, kpc. 

2-7: U, (Ty, V, ay, W, aw — Velocities transformed to LSR frame from the measured proper motion and 
errors. At each distance, od is taken as 0. E.g., At 10 kpc, a\j, ay, aw, are 8, 8, 4 km s _1 , from errors on 
Mra, A^dec, and V Hc iio only. 

8,9: UC, VC = Velocities (LSR), as expected from circular rotation at 220 km s _1 . WC = 0. 
10: VRh = Radial vel (Heliocentric), expected from circular rotation. 

11: VTh = Transverse vel (Heliocentric) across line of sight, expected from circular rotation. 
12: VRq = Measured radial vel (Galactocentric). km s -1 expected for circular rotation. 
13: VCg = Measured circular vel(Galactocentric). 220 km s _1 expected for circular rotation. 
14: V P = Peculiar velocity = [(U-UC) 2 + (V-VC) 2 + W 2 ] ' 5 = [VR G + (VC G - 220) 2 + W 2 ] - 5 
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Fig. 1. — Astrometric Position of the Black Hole over a timespan of almost 10 years. 



Velocity Components 




Fig. 2. — Measured velocity components (solid lines) vs. distance, compared with expected 
values for an object in circular rotation about the Galactic Center (dashed lines). 
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Fig. 3. — Total peculiar velocity of GRS 1915+105 plotted at various assumed distances. 
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Fig. 4. — Galactocentric Orbit of GRS 1915+105 for the past 1 Gyr, (dashed blue line), with 
the most recent 250 Myr shown as the solid red line. 



